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The electromagnetic response of graphene in a magnetic field is studied, with particular emphasis 
on the quantum features of its ground state (vacuum). The graphene vacuum, unlike in conventional 
quantum Hall systems, is a dielectric medium and carries an appreciable amount of electric and 
magnetic susceptibilities. The dielectric effect grows rapidly with increasing filling factor v in such 
a way that reflects the 'relativistic' Landau-level characteristics of graphene as well as its valley and 
spin degeneracy. A close look into the dielectric function also reveals that the Coulomb interaction 
is efficiently screened on the scale of the magnetic length, leading to a prominent reduction of the 
exciton spectra in graphene. In addition, an effective gauge theory of graphene is constructed out 
of the response. It is pointed out thereby that the electric susceptibility is generally expressed as a 
ratio of the Hall conductance to the Landau gap. 
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I. INTRODUCTION 

There has recently been increasing interest, both 
experimentallji^i^ and theoreticallj«^i^iii^i2iiiS, in a 
'relativistic' condensed-matter system, graphene which 
is a monolayer of graphite. Graphene is a gapless pla- 
nar semiconductor, in which low-energy electronic trans- 
port is essentially governed by massless Dirac fermions 
with effective speed of light wp ~ 10^ m/s « c/300, and 
thus provides a special opportunity to study relativistic 
quantum dynamics in condensed-matter systems. Ex- 
periments have revealed a number of exotic transport 
properties of graphene, such as an unusual sequence^'^ 
of magneto-oscillations and of the quantum Hall (QH) 
effect, that are characteristic of Dirac fermions. 

It has long been known that Dirac electrons in 2+1 
dimensions lead to peculiar quantum phenomena of frac- 
tional fermion numberJi: Fractional charge and current 
of abnormal parity, as summarized by the Chern-Simons 
(CS) term, are induced in the vacuum in response to an 
applied field. The induced CS term, or a parity anomaly, 
is associated with a spectral asymmetry and an index 
of the Dirac Hamiltonian, and is also tied to the chiral 
anomaly in 1-1-1 dimensional^. The one-half degeneracy 
of the lowest Landau level observedi*^ in the graphene 
QH effect is a manifestation of fermion number fraction- 
alization, although complete cancellation of the net vac- 
uum charge and current takes place in graphene (which 
involves a parity pair of fermions). 

The purpose of this paper is to explore possible sig- 
natures of relativistic quantum field theory in the low- 
energy physics of graphene. We study the electromag- 
netic response of graphene in a magnetic field at integer 
filling factor z^, with particular emphasis on the quan- 
tum features of the ground state with v — Q, the vacuum 
state. For graphene, unlike standard QH systems, the 
vacuum state is a dielectric medium and carries an ap- 
preciable amount of both electric and magnetic suscep- 



tibilities (aejQfm) over the entire range of wavelengths; 
this reflects the presence of the 'Dirac sea'. The dielec- 
tric effect grows rapidly with increasing filling factor v 
in such a way that reflects the 'relativistic' Landau-level 
characteristics of graphene as well as its valley and spin 
degeneracy. A close look into the dielectric function re- 
veals that the Coulomb interaction is efficiently screened 
on the scale of the magnetic length, leading to a promi- 
nent reduction of the exciton spectra in graphene. 

In addition, we construct, out of the response and via 
functional bosonization, a low-energy effective gauge the- 
ory of graphene in a magnetic field. It is pointed out 
thereby that the electric susceptibility ac is generally ex- 
pressed as a ratio of the Hall conductance to the Landau 
gap. 

In Sec. II we study the electromagnetic response of 
graphene in a magnetic field, with emphasis on projec- 
tion to Landau levels and regularization. In Sec. Ill we 
construct an effective theory of graphene. In Sec. IV we 
discuss the effect of polarization by treating the Coulomb 
interaction in the random-phase approximation (RPA). 
Section V is devoted to a summary and discussion. 



II. ELECTROMAGNETIC RESPONSE 

Graphene has a honeycomb lattice which consists of 
two triangle sublattices of carbon atoms with one electron 
per site. The electrons are described by a two-component 
spinor [U, VY with fields U and V residing on each sub- 
lattice. In the tight-binding approximation with nearest 
neighbor interactions the spectrum of electrons becomes 
linear at the two inequivalent Fermi points {K and K') 
on the corners of the Brillouin zone. In the continuum 
limit the system is described by a low-energy effective 
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Hamiltonian of the for: 
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H 



'H± = i;f (crilli + (72112 ± TOCTa) — e^o, (2-1) 

with couphng to external electromagnetic potentials 
— (AijAo) introduced through 11^ = —idi + eAi 
[i = (1,2) or (x,2/)]; v-p ~ 10^ m/s is the Fermi ve- 
locity. Here ijj and x stand for the electron fields near 
the K and K' points (at wave vectors =Fq) and have 
the structure ipa = (V'1,'02)' oc (C/k-q, Vk-q)' and Xa = 
(Xi,X2)' oc (-Vk+q,C^k+q)*, where C/k = J7k(t), e.g., de- 
notes the Fourier transform of ?7(x, t) for short. 

For generality we have introduced a tiny 'mass' gap to, 
which works to lift the degeneracy of the lowest (n = 0) 
Landau level alone. The possibility of such a mass gap 
has been discussedZ.'-^-'?'^" in connection with the observed 
lifting^ of the degeneracy of the n = 0, 1 levels in high 
magnetic fields. Nonzero mass m ^ spoils the pseu- 
dospin (or sublattice) SU(2) symmetry of i7, which ro- 
tates {tp, x)- Actually, we keep to to reveal the parti- 
cle/hole character of the n=0 levels, as remarked below, 
and set to ^ eventually. 

For clarity, we suppress the electron spin, which is 
treated as a global SU(2) symmetry of H by doubling the 
fields, and x° with a = (t,i). The SU(2)-breaking 
Zeeman coupling is considerably weak, compared with 
the Coulomb interaction, for graphene and is ignored in 
the following. 

The Coulomb interaction is written as 



H 



Coul 



2 ^ 
p 



(2.2) 



where pp is the Fourier transform of the electron number 
density p = ip'^tp -\- x^X: "^p — 27ra/(eb|p|) is the Fourier 
transform of the Coulomb potential w(x) = a/(eb|x|) 
with the fine-structure constant a = e^/(47reo) ~ 1/137 
and the substrate dielectric constant eb. In this pa- 
per we shall focus on the long-wavelength properties of 
graphene and consider H + TJ^oui ^^j. j-^g^gj^. Hamil- 
tonian; other short-distance corrections^, including the 
oscillating piece oc e^^^'i-^ of the Coulomb interaction, 
will thus be ignored. 

Let us place graphene in a strong magnetic field and 
first study how the electrons in graphene respond to 
weak external potentials Afj_(x). To this end we set 
Ai{x) (x)-l-Ai(x, t) in the Hamiltonian ([2T|) . where 

the vector potential = B{—y,0) supplies a uniform 
magnetic field Bz = B > normal to the graphene plane. 
We turn off ijCoui ^^j, ^j^g moment. 

When Af^ — 0, the eigenmodes of H are Landau levels 
of ■0 and X of energy 



(2.3) 



labeled by integers n = 0, ±1, ±2, . . ., and Px (or yg = 
i'^Px with the magnetic length i = 1/ VeB) ; ujc = vp jl 



is the basic cyclotron frequency. Here s„ = sign{ri}(= 
±1) specifies the sign of the energy e„. 

For each n ^ ijj and x have the same spectra and 
the positive- and negative-energy (i.e., n > and n < 0) 
levels are symmetric in structure. [This is a consequence 
of the conjugation symmetry 7i+ = —aj,{TL-) with 
Aq ~ 0.] The n = spectra depend on the sign of to. 
Let us take to > 0. Then the n = level of il) has negative 
energy eo_ — — fpTO while that of x has positive energy 
eo^. = vpm. These n — levels represent holes and 
electrons via quantization; and this hole/particle char- 
acterization persists even in the limit to ^ where the 
n = 0::p levels become degenerate. These n = 0± eigen- 
modes have components only on each separate sublattice; 
see Appendix A for details of the wave functions. With 
the electron spin taken into account, each Landau level 
is thus four-fold degenerate, except for the n — Q± levels 
which are doubly-degenerate. 

To make this Landau- level structure explicit, it is 
useful to pass to the |n,?/o) basis, with the expansion 
V'(x,t) = Y.n,yo^^\'^^yo)'^n{yo,t). (From now on, we 
shall only display the ip sector since the x sector is ob- 
tained by reversing the sign of to.) The translation is 
simpl e ^"i"— '! In the |n,yo) representation the coordinate 
x = {xi^X2) of ^^(x, t) is split into the center coor- 
dinate r = {ri,r2) — {i£'^d/dyo,yo) with uncertainty 
[^i,''2] = and the relative coordinate X = (Xi,X2) 
with [Xi,X2] = -iP which are matrices in the Landau- 
level index. The Hamiltonian H thereby is rewritten as 

H ^ Jdyo J2 ^n(2/o,i){ + eVnn'}lp7i'{yo,t), 

n,n' — — oo 

Vnn' = VpAmi' + f F ~ (A)rm'- (2.4) 

Here the electromagnetic coupling Vnn' is a matrix in the 
level index and a function of r = {ii'^d/dyo, yo): The field 

(A)«„'(r,t) = ^e^P-"-w(p)e-^''p' (Ao)p (2.5) 
p 

is expressed in terms of the Fourier transform {Ao)p{t) 
of Aq(x, t) and the coefficient matrix 



Onn' (p) 



C+ C+ /|„|_l>'|_l(p) 
^SnSn'C" C~, /|„l,ln'l(p) 



(2.6) 



with = y^l ± VYtn/tn, where 
/fcn(p) = 



n\ /i£p\k-n 



i^V) (2.7) 



for k > n; p = py+ ip^', for k < n one may set k ^ n 
and p Py— ipx- Similarly, Ann' is defined by Eq. (|2.5[) 
with (^o)p ^ {-^x + iAy)-p and 

5rm'(p) jnn'(p) = 5 C+ C^' - 1, |n' | (p) ; (2.8) 
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■^in' is defined with {Ax-iAy)p and j^„,(p) = j„'„(-p). 
It is clear from Eq. (|2.4p that (Ao)p is coupled to the 
charge density P-p{t) = J (Pxe^^'^ ip'^ip rewritten as 



5rm'(p) / rfyO^ie'P-'-Vn'; (2.9) 



n,n' 

analogously for the current J{x) = ^ '4'H,<^i + ^(^2)'^ cou- 
pled to A ^ Ax + lAy. 

Even a weak potential (x) causes mixing of the Lan- 
dau levels, and its effect is calculated by diagonalizing H 
with respect to the true levels {n} by a suitable U(oo) 
transformation ip^{yo,t) = J^m'^rimir^t) iprniyo,t). 
This yields H = J dyoEni^'^Vn'^nn^n with H« = 
G{H - idt)G-^ diagonal in the level index. To 0[Al), 



with 



^r?^ = ^El{^»fe ^TtT ^'-^ - (2-10) 

where i9t acts on Vkn- This O(y^) term embodies the 
{pp) and {JJ'^) correlation functions and constitutes the 
linear response of the system. Actually it only takes 
account of inter-Landau-level transitions due to the cy- 
clotron modes; they govern the response of graphene at 
integer filling factors, on which we focus throughout this 
paper. For graphene at noninteger filling one encounters 
intra-level (collective) excitations,^^ which are less sen- 
sitive to long- wavelength probes and are determined by 
diagonalizing the Coulomb Hamiltonian ijCoui projected 
to a given Landau level. 

Let us now consider the long-wavelength part of the 
electromagnetic response, which, in view of gauge invari- 
ance, takes the form 



-Pn \ e^Ekf + Vn \e\A^2f + ■■■, (2.11) 

where Ek = -<9feAo - dtA^ and A12 = diA2 - 92^1; 
gOi2 _ 2. Actual calculations are simplified if we take 
A^(x, t) uniform in xi, i.e., A^ljj^t). Then the CS term 



dn{yo,t) = (V''^)!L„(yo, *) with n > to denote electrons 
and holes over the vacuum. The Hamiltonian then reads 

/oo 00 
dyo[j2hn aia„-J2 h-nAdn+ a] , (2.13) 

n=l ri=0 

00 

A = (5o ^ /i-„, (2.14) 

n=0 

where = <5(yo = 0) = Lx/i^-Ki"^) stands for the degen- 
eracy of each Landau level; = J dxi. 

Here A represents the quantum response of the vac- 
uum. Apparently it is indeterminate because of the sums 
over an infinite number of Landau levels. Note, e.g., that 
7„ = e^^^ > for all n. This would naively mean that 
the vacuum, the Dirac sea, carries an infinitely large Hall 
conductance, which is physically unacceptable. To ob- 
tain a sensible answer one has to define the sum care- 
fully. For regularization let us truncate the spectrum 
to a finite number (2A^ + 1) of Landau levels {n} with 
—Cat < En < Cat, and let ^ 00 at the endJ^ 

It is instructive to see why and how this regularization 

/i2\ 

works physically. Let us note that in Vn of Eq. (|2.10[) 
the virtual (n k n) transition and the related 
(k ^ n ^ k) transition contribute equally but in op- 
posite sign. For definiteness, we denote by F„ the 0{V^) 
response, 7n, '7n, . . ., collectively, and write it as the 
regularized sum F„ = J2k=-N o'^^'" ^be contribution 
from the n ^ k ^ n subprocesses. Then, in the 
sum Fn cancellation takes place among a majority of 
terms, owing to the antisymmetry 



(2.15) 



which implies pair-wise cancellation between the n — > 
k ^ n and k n ^ k subprocesses. Fermi statistics is 
thus naturally taken care of in the regularized sum. Such 
a pair-wise cancellation of virtual processes is a basic 
property of Berry's phasei^ii^; and it is also an exact 
property^*^ of the density and current response functions 
(with also the Coulomb interaction included), as is clear 
from their spectral representations. 

The vacuum response is now written as a regularized 



e'^'^^A^d^Ax ^ {AyAx - AxAy), the Ef - [A^f + (4)^ sum A = <5o E^Lo '^-^ P 



term and the (^12) ~ {dyAx) term are derived from 
the 0{dt), 0{d^) and O(p^) parts of the {JJ'^) response 
function, respectively. A direct calculation yields 



that the regularized sum of F„ over all the levels vanishes, 



cast i^™'^ in another suggestive form 



0. This fact allows one to 



2wf 



nf 

-N 



n— 



nf N 
n= — N n>nf 



(2.16) 



(2.12) 



Note that 7„ is independent of n while /3_„ = — /3„ (with 
/3o^ = =f|w|£''/wf) and = (with 770 = 0). 

We are now ready to introduce the quantum vac- 
uum as the Dirac sea with all negative-energy electron 
states occupied. Let us rename a„(?/o,i) = V'n (j/Oj^) and 



Here we have written the response in a slightly general- 
ized form: The left-hand side is the response of a many- 
body state with Landau levels occupied up to rt = nf ; the 
choice nf = 0_ thus yields F™'^. This formula expresses 
the quantum response in terms of an asymmetry in the 
spectrum of the occupied and empty levels, weighted with 
the response Fn per level. 
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Via regularization the calculation of the response is 
modified. A close look into the matrix elements in 
Eq. (|2.8p shows that /3„ and 7„ come only from virtual 
transitions to the adjacent levels {n ^ n±l) and the re- 
lated ones across the Dirac sea {n — {n ± 1}). As 
a result, the expressions for 7„ and in Eq. (|2.12|) 
are valid for |n| < N — 1 while the boundary contri- 
butions 7-Ar and /3_jv have to be calculated from the 
—N ±(A^— 1) ~N processes. For the result is 



the Dirac sea. This reproduces G(0) as the sum 



JV-l 



6G(0) = l+J2i^nTl-V^f 



(2.23) 



fe=i 



(2.17) 



which then makes A""" finite, 



AT 



2^1-n^^eiSo = L^^^^ 

n=0 



7Z^- (2-18) 



This implies that the vacuum would have Hall conduc- 
tance equal to half of a filled level, CTnaii = ^e^/ft., and 
that the Hall effect would arise without real electrons or 
holes. This nonzero A''' is attributed to a spectral asym- 
metry due to the n = 0_ level if one uses Eq. (|2.16p (since 
j-N — JN — 0), as normally implied by an index theo- 

1 9 

ren>». 



Formula (|2.16p also leads to the same result (|2.20p . 
Note that this formula admits a further generalization: 
Instead of the symmetric spectral cutoff one can equally 
well truncate the positive- and negative-energy spectra 
independently by choosing —TV < n < TV' so that there 
are {N -I- A^' -I- 1) Landau levels in total. Clearly f3-N and 
7_Ar are unchanged if one takes N' > A^. 

With this generalization in mind let us now exam- 
ine the magnetic susceptibility of the vacuum, a™'^ = 

{J2n=i''l-n + rjO-) ^o/ Lx- A look into the matrix ele- 
ments again reveals that the expression (|2.12p for r/n is 
valid only for \n\ < N — 2. Also a direct calculation 
shows that the virtual transitions from the bottom of the 
Dirac sea to positive-energy states combineS^ to vanish 
(for A^ ^ oo) once one takes A^' > A^ -I- 2, and the re- 
sponse is independent of the way one sends (A^, A^') oo. 
We therefore set A^' > A" -I- 2 to calculate the boundary 



The physical mechanism underlying the vacuum term Ty.jv + 77_(Ar_i) and obtain 



Hall effect was discussed earlier-iSi^. Actually, the x sec- 
tor contributes to A'*' in opposite sign and the vacuum 
Hall effect does not take place for graphene. 

In contrast, the vacuum effect survives for suscepti- 
bilities Pn and r]n that change sign with n. For /?„ the 
boundary contribution again increases rapidly with A^, 



3A| 
47r Wc 



(2.24) 



with 



/3^jv« V2(^V«f) 



±(A^- 



iA) 



3/2. 



(2.19) 



where A — ^m^£^, and makes the vacuum electric sus- 
ceptibility a™'^ 



^ n + A V 

n— 1 

i( 



n + \ + 



A/3 
VnTA 



5 ^ i A 



3/2 



iY.n=i + /^o_ ) 5q/Lx finite, 
3e2 



with 



2tt uj. 



N-l 



{G(A)-ix/A}, 



-C(-i) + tC(i) + 

0.2079- A -0.2434 4 



(2.25) 



(2.20) 



G(A).-E{^/^-i§,} 



-C(-i)-AiC(i) + 

0.2079 -I- A • 0.2434 



^/n + X■ 

3/2 



Let us next consider the Aq response in /i„, which 
probes the charge of the system. The observable charge 
is a difference between the B and B = cases. Such 
an induced charge is detected by a normal-ordered charge 
operator : 'ip^i'ipn ■ — \ VPn^ V'n] and a symmetric spectral 
cutoff, with the resulti^ 



AQ"'^'= = - 



2 2nP 



e2B/(47r) 



(2.26) 



(2.21) 



Here we have used 

N-l 



'n - 



{N-\f'^^ _0.2079, 



E 

n=l 

JV-1 

Y^l/^^2^N -\ = C(5) « -1-4604. (2.22) 

n=l 

It is instructive to derive this a™'^ by selectively sum- 
ming up the virtual (— n ^ n ± 1 — > —n) processes from 



per unit area. Alternatively, this induced charge is de- 
rived from A"*" in Eq. (|2.18p if one notes that the variation 
of AQ under the change B ^ B + A12 is read from the 
CS coupling. It is thus clear that there is no induced 
vacuum charge in graphene. 

Similarly, one can calculate the observable vacuum en- 
ergy density evac = ^rac— ^vac° from the e„ terms in hn via 
proper regularization; see Appendix B for details. Here 
we only remark that this Cvac provides another way to 
derive a™'^ : Setting B B + A12 in evac and expanding 
it to second order in A12 yields precisely the regularized 
expression (|2.24p for a^^. 
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We are now ready to write down the electromagnetic 
response of graphene. Let us suppose that the electrons 
fill up an integral number v of Landau levels, with uni- 
form density {p) = p = z//(27r£^); the case of holes is 
treated analogously. We take both spin and pseudospin 
into account, and write v — in terms of the filling 

factors Vn of the nth level [0 < < 4 for n > 1 and 
< < 2 ] . The long- wavelength response of graphene 
is then summarized by the Lagrangian 

La = peA^~e'ep]^e^'''P A^d.Ap 

+ iaeME2-ia„,(i.)(Ai2)2 (2.27) 



with the susceptibilities 



ctc{v) 



3e^ 
27rWc 



oc 

[^i/„V^ + 4G(0) 



{vl/2)aM- 



(2.28) 
(2.29) 



For clarity we shall from now on display formulas with 
m — > 0. The vacuum susceptibility a™'^ is almost com- 
parable to the contribution of a single filled n — 1 
level, i.e., 4 x 0.21 vs 1. Curiously, in contrast, the 
n — Q± levels hardly contribute to ac{v)\ their tiny frac- 
tion cx (1^0+ ~ 2)(V2/3) \m\t in ac shows that a mass gap 
works to reduce Ue, as is clear intuitively. 



III. LONG- WAVELENGTH EFFECTIVE 
THEORY 

hi this section we derive an effective gauge theory 
of graphene in a magnetic field. It is constructed, via 
functional bosonization^-, so as to reproduce the re- 
sponse (|2.27p . Such an effective theory has some ad- 
vantages2^i2i: It allows one to handle the Coulomb in- 
teraction exactly. It readily admits inclusion of vortices, 
and is thus applicable to the description of the fractional 
QH effect as well. It furthermore unifies^^ the composite- 
boson^*' and composite-fermion^^ descriptions of the frac- 
tional QH effect. 

Applying-2S, the standard procedure of functional 
bosonization to the response (|2.27p yields an effective the- 
ory of a vector field 5^ ~ (607 ^2), with the Lagrangian 
to 0(^2), 

Lb - T^b^e^^'^^d^bx + ^ho 



2Pp^ 



and the effective cyclotron frequency 



3 T.n=0^nV^+^G{Qy 



(3.2) 



where b^i, = d^b^-d^b^^; Pp = v/(2tt) and v = Y.°^=q ^n- 
The Coulomb interaction and coupling to are incor- 
porated into this theory, with the Lagrangian 

Leff[6] - ~eA^ e^^'^dM - \ 6bi2 v Sh^ + Lb, (3.3) 

where Sbi2v5bi2 = J (Px'Sbi2{x) v{x—x')6bi2{x') for 
short and dbi2 — ^12 — P- 

This effective theory reproduces the original re- 
sponse (|2.27l) . and actually more: The AedA, and 
(^12)^ terms acquire a common kernel so that, e.g., 
E2 ^ E /C E, with 

JC = c^e'ff/K'ff + df - {ujosfpvp + anja,) V^]. (3.4) 

This shows that the Coulomb interaction Vp ^ e^/jpl 
substantially modifies the dispersion of the cyclotron 
mode at long wavelengths p ^ 0, 



1 V 
W(p) W CJoff + - {fpVp + 

Z ZLOcS 



(3.5) 



For graphene the Landau levels are not equally spaced 
and the excitation gaps depend on the level index n or 
V. When the nth levels are occupied, i.e., at = 4n + 2, 
the minimum activation gap is Awc"'' = Cn+i — (n = 
[s/n -|- 1 — ^/n) LOc] numerically, Awc*^^ = l.Oo^c, Ao^c"""^ ~ 
0.4142 Wc, Aw^^' w 0.3178 Wc, Aw^^^ « 0.2679 Wc, for v = 



2, 6, 10, 14, respectively. The Wcff in Eq. (|3.2p represents 
such an excitation gap. 

For 1/ = 2, 5(2) = 1/{6G(0)} « 0.802, which comes 
solely from the vacuum fluctuations, appreciably devi- 
ates from the number 1.0 quoted above. For v = Q, 
g{&) ~ 0.4141; for u = 10,14,..., g{u) almost precisely 
agrees with the quoted numbers, and the agreement be- 
comes exact for — > 00. It is somewhat surprising that 
an effective theory constructed from the long-wavelength 
response alone gives an excellent description of the exci- 
tation spectrum. 

In other words, the effective theory reveals a general 
link: The susceptibility Uc of a QH system is essentially 
determined by a Hall conductance and the Landau gap. 



/Ac.("). 



(3.6) 



This implies that q;o(j^) rises linearly with v for typical 
QH systems with equally-spaced Landau levels while it 
grows more rapidly for graphene. 



(3.7) 



Actually, formula (|3.6p is a condensed-matter realization 
of the related fact'^® in relativistic field theory that in 
2-1-1 dimensions the gauge field acquires a mass from the 
CS term in the presence of the Maxwell term ^ E^ . 



IV. POLARIZATION FUNCTION 

We have seen that the vacuum component a™'^ consti- 
tutes an appreciable portion of the susceptibility aci^v). 
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In this section we examine the effect of polarization at 
shorter wavelengths and extract some observable conse- 
quences. We shall see that the screening properties of 
graphene in a magnetic field are substantially different 
from those of graphene without the magnetic field, dis- 
cussed recently^-. 

Let us turn off the Coulomb interaction for the mo- 
ment and calculate the polarization function P{p, lu) = 
—i{p{x)p{y))^'^', where F.T. stands for the Fourier 
transform. To this end we denote by P„(p,w) = 
—i {n\Tp{x)p{y)\n)^''^' the polarization function for the 
filled nth Landau level \n) — {|n,?/o)} and write it as a 
regularized sum P„(p,ct;) = J2k=-N ^niP^^) O'^sr com- 
ponents coming from the virtual (n — > fc — *■ n) tran- 
sitions. As before, we truncate the spectrum to a finite 
interval -N <n,k<N. From Eq. (pTTO)) one finds 

F„^p,c.) = -{^— + ^—\a';,{p), (4.1) 



27r^2 



e 2 ■ 



5nfc(p)gfc«(-p), (4.2) 



where tkn = ^k — ^n', "'^(p) is the spectral weight. 

One can calculate the response P{p, w) = Pn{p, ^) 
by selectively summing P^ over occupied levels {n} and 
unoccupied levels {k} for a given graphene state. In prac- 
tice, it is advantageous, especially for numerical evalua- 
tions, to calculate P„ = SfcL-Ar Pn by summing over all 
levels {k} since Fermi statistics is automatically taken 
care of in the regularized sum, as noted in the previous 
section. We focus on the real part of P„(p,tt') in the 
static limit uj ^ below. For the ip sector the result is 



Pn{p,0) 

Up') 



3p2 



2tt cu. 



S„Cn(P^), 



e 2 



(4.3) 



|n|e(i^V,|n|), (4.4) 



with 



E(n-l)! x" ^ 
(n + fc-1)! k 



k=l 

-E 

k=l 



n+k 



ik)\2 



(n — fc)! X 



I ^k-l r 



71! 

k 



k 



(jik) ^ 



V n 



(fe) 



k ~^ ^n-k- 



n — k 
c"^i 1 



n! n Anx 



(i„ - L„_i)2 (for n > 1), (4.5) 



and f(x,n = 0) = 0, where lI'^-' — L^\x) for short; 
f(x,n) are normalized so that ^(0,n) = 1 for 1 < n < 
N ~ 1, with their reference to the cutoff N suppressed. 
Note that P_„ ~ —Pn\ in particular, Pq = 0, i.e., the 
n = 0± levels do not contribute to polarization (for m 
0). 



One can now write P(p, 0) or the static susceptibility 
function ac[p] = — (e^/p^) P(p, 0) for graphene as 

«c[p] = \ Y.^nUp')+^r'{p% (4.6) 

n>0 

where C™=(p2) = - E^iln(p^). This a^[p\ is correctly 
reduced to in Eq. (|2.28p for p ^ 0. 

Let us now turn on the Coulomb interaction Vp = 
27rQ!/(eb|p|) and study its eflFects in the RPA. The RPA 
dielectric function is written as 



e(P,t^) = ^-VpP{p,uj), 



(4.7) 



and the RPA polarization function PF;pA(p,ti^) = 
P(p, w)/e(p, w) is screened through e(p,w). The static 
dielectric function then reads 



e(p,0) = l 



3a £\p\ 



CbWF %/2 



n>0 



where a — e^ji^Tito) w 1/137 and wp 



(4.8) 

c/300. With 



air on one side of the graphene plane and Si02 on the 
other, the unscreened dielectric constant is estimated^ to 
be Eb ~ 2/ (e^^J + egjQj « 1.6 with eah- « 1 and esiOs ~ 4; 
this yields an estimate of the coefficient 



3a/(ebWF) — iV2a/ {tx^iojc 



4.1 



(4.9) 



In Fig. 1 we plot a^lp] and e(p,0) — 1 for v ~ 
2, 6, 10, 14. Note first that there is no screening at 
long distances, e(p, 0) ^ 1 for p ^ 0, as is typical of 
two-dimensional systems. However, e(p, 0) grows rapidly 
with IpI and becomes sizable for |p|^ ~ 1, with its peak 
value e(p, 0) w 3.4 ^ 20 for v — Q ^ lA. For comparison, 
analogous plots for a standard (GaAs) QH system2£ (per 
spin) are shown in Fig. 2. Both ae[p] and e(p,0) — 1 
look qualitatively similar in the two cases, but there are 
clear differences: First, for graphene even the vacuum 
(with V — Q) has a finite polarization, which decreases 
only gradually with increasing |p|^; this reflects the fact 
that the Dirac sea has quantum fluctuations of a whole 
range of wavelengths. Second, for graphene ac[p] and 
e(p,0) grow fast with v. This reflects the Landau-level 
characteristic of graphene; the quantum effect grows as 
the level gaps become narrower. 

To compare graphene and GaAs systems let us note the 
following: For GaAs Wc = eB/m* « 20P[T] K is con- 
siderably smaller than for graphene while the substrate 
dielectric constant eb ~ 12.9 3> (eb)graphcnc ~ 1-6 so that 



(Cb '^c)GaAs/ (cb w 0.4 VP[T]. (4.10) 

Figure 2 indicates that the effect of polarization is gen- 
erally not quite appreciable for the GaAs system, with 
e(p,0) « 1.7 ~ 2.2 at |p|^ - 1 for i/ = 1 - 3 ff we 
take 3\/2 a/(eb^'^c) ~ 4, which, in view of Eqs. (|4.9p 
and (j4.10p . would overestimate e(p,0) numerically over 
a typical range of the applied magnetic held _B > 10 T. 
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FIG. 1: (a) Susceptibility Qe[pl for graphene at integer fill- 
ings u, plotted in units of /{2iyloc); n refers to the highest 
occupied Landau level, (b) Dielectric function e(p, 0) — 1 is 
plotted in units of 3a/(eb vf) ~ 4.1 so that the peak value of 
e(p, 0) ranges from 3.4 to 20 for = ~ 14. 



e(p,0)-l [3V2a/ebf^J (q^^s) 
.3[ 




2 4 6 8 



FIG. 2: Susceptibility [p] and e(p, 0) — 1 for a GaAs system 
(per spin). The peak value of e(p, 0) ranges from 1.7 ~ 2.2 
for I' = 1 ~ 3 if one takes 3\/2 Q!/(eb^'^c) ~ 4. 



One can derive^ the excitation spectrum corrected 
by the Coulomb interaction from the RPA response 
^'rpa(p,w) = P{p,uj)/e{p,uj). Let us isolate from 
P(p, uj) one of its poles at a; ~ ekn — — £n and set 
e{p,uj) — 1 — Vp P{p, uj) 0. This fixes the pole position 
of the RPA response function, e^^f" = Cfc — + ^^k,n{p) 
with 

Aefe,„(p) « ^z.g(p)€|p|/{2e(p,0)}, (4.11) 

^g(p) =^9nk{-p)gkn{p)e-''/x, (4.12) 

where x = ^£^p^; the sum is taken over possible (/c,n) 
that lead to the gap — e„. We have parametrized 
Aefc.„(p) so that i'g{0) corresponds to z/ = in the 

excitation spectrum (|3.5p of the effective theory. Note 
that Aefc^„(p ^0) = so that the excitation gap 
Efc — remains unshifted for p — > 0. The excitation 
gap Awc"^ = f^c: in particular, is realized via the 0± — > 1 
or — 1 ^ 0+ transitions, yielding 

t'g(p) = 2e"*^'p' for v = Q,2. (4.13) 

The next gap AtJc^-* QAAlo^ aX v = Q is due the 1^2 
transition, giving fg(p) = (1 + \/2 — xl^/2)'^e^^ with 

X = \ £'^p^. Also the gaps Aw^"^ = {^n + 1 - -/n)wc 
open up at — 4n + 2 via the n ^ n + 1 transitions with 
i'g(O) = (v^ + VfT+T)^. Actually, these numbers fg(0) 
and 1/ = 4n + 2 agree within 3 percent for n = 1 and 
the agreement improves rapidly for larger n. This shows 
that the excitation spectrum (|3.5p in the effective theory 
is practically in good agreement with the RPA result. 

In Fig. 3 (a) and (b) we plot the exciton spectra 
Aefc^„(p) near the energy gaps Wc, 0.41 cjc and 0.32 tJc. 
The unscreened spectra (dashed curves) are generally 
peaked around |p|^ ~ 1 and are efhciently reduced to 
the screened spectra (solid curves) through e(p, 0) also 
peaked at \p\l ^ 1. Near the gap o^c at = 0, 2 the 
spectrum is about one third of the unscreened one in 
magnitude. The effect of screening is far more promi- 
nent at higher v. Near the gap 0.41 lu^ dX v = & the 
spectrum is about one order of magnitude smaller than 
the unscreened one. Near the gap 0.32 Wc at = 10 the 
spectrum is even lower than the v = & spectrum although 
the unscreened spectra are opposite in magnitude. 

For comparison Fig. 3 (c) shows the exciton spectra in 
the GaAs QH system'^*' near the energy gap ojc — eB/m* 
at V = 1,2. There the spectra are reduced roughly by 
50 percent or less. Clearly the effect of screening is not 
quite important for GaAs systems. 

V. SUMMARY AND DISCUSSION 

In this paper we have examined the electromagnetic 
characteristics of graphene in a magnetic field to explore 
possible signatures of quantum field theory in the low- 
energy dynamics of graphene. The key features that 
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!/ = 10 (n = 2) 
i/ = Q (n=l) 



IpI^ 
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FIG. 3: (a) Excitation spectrum near the gap ujc at v = 
and 2, plotted in units of a/ebi- It is reduced from the 
bare spectrum (dashed curve) via screening, (b) Excitation 
spectra near the gap 0.41 Wc and 0.32 cOc at v = G and f = 10, 
respectively. The unscreened spectra (dashed curves) are one 
order of magnitude larger than those depicted in the figure, 
(c) Excitation spectra for a GaAs system per spin, with the 
choice 3\/2a/{e-h£ujc) ~ 4. 



distinguish graphene from standard QH systems are (1) 
the quantum nature of the graphene ground state (vac- 
uum) and (2) the relativistic pattern of Landau levels 
(oc ±-\/|n[) with gaps decreasing with the level index n. 

The graphene vacuum is a dielectric medium with elec- 
tron and hole pairs created (from the depth of the Dirac 
sea) in response to an applied field. In particular, the 
electric and magnetic susceptibilities (aejCtm) and the 
dielectric function e(p, uj) — l are nonzero for the vacuum 
over a whole range of wavelengths and grow prominently 
with increasing filling factor in a way that reflects the 
Landau-level characteristics (2) of graphene. 

As for experimental verification, detection of the inter- 
Landau-level cyclotron-mode excitations, or excitons, 
from the v ^ ground state through light absorption or 
reflection would be a direct signal for the quantum nature 
of the vacuum. A detailed study of the exciton spectra 
from the states with = 0, 2, 6, . . . would also reveal the 
effect of screening around |p|£ '--^ 1, which is expected to 
be sizable for graphene, as discussed in Sec. IV. 

In Sec. II special emphasis has been placed on the need 
for regularization in calculating quantum corrections for 



graphene. We have pointed out that the regularized (lin- 
ear) response is written as an asymmetry in the spectra 
of the occupied and unoccupied states, weighted with 
the response per state; see Eq. (I2.16p . For the induced 
CS coupling this spectral asymmetry is related to the 
?7-invariant of the Dirac Hamiltonian (i.e., the difference 
in number between positive and negative eigenvaluea^^) , 
although it vanishes for graphene eventually. For the 
susceptibilities or the polarization function, in contrast, 
there is no such underlying index theorem but the spec- 
tral asymmetry yields nonzero results for them, with the 
contribution from each level adding up. In this connec- 
tion, some peculiarities of the n — Q± Landau levels 
would be worth noting. While they carry the (normal) 
Hall conductance zte'^/h per level, they hardly contribute 
to susceptibilities and (for m — > 0). As a result, 
it is essentially the vacuum state that would govern the 
dielectric property of graphene for < 2; the effect of 
polarization would be insensitive to carrier densities over 
this range; see, e.g., Eq, (|4.13p . 

In Sec. Ill we have seen that an effective gauge the- 
ory of graphene, constructed from the long-wavelength 
response, gives an excellent description of the excitation 
spectrum. In this effective theory graphene and other 
QH systems look quite similar, except that the effective 
Landau gap Awc ~ Weff in Eq. (|3.2p now depends on v 
for graphene. We have thereby noted a general relation 
etc ~ fHaii/ij^eff which relates the susceptibility ae to the 
Hall conductance and the Landau gap. 
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APPENDIX A: FREE DIRAC ELECTRONS 

In this appendix we summarize the eigenmodes of the 
Dirac Hamiltonian in a magnetic field. Let us focus on 
the ip field governed by ?{+ = vp (ctiHi -|- a2ii2 + rna^) 
in Eq. (HH) with A, = B{-y,0). The spectrum 

of ijj consists of a n infinite towe r of Landau levels of en- 
ergy e„ = Sn Wc \/\n\ + wP-P- jl labeled by integers n = 
0,±1.± 2, • • •, with I = 1/%/eB and s„ = sign{n} = ±1. 
The positive-energy {n > 0) and negative-energy (n < 0) 
eigenmodes are neatly cast in the unified form 



i^nyo (x) 



1 

V2 



C.^ 0|nhl(2/-yo) 
SnC~ (t)\n\{y - yo) 



1 



V2^ 



(Al) 

where j/o = ^'^Px and = a/I ± vprn/en, in terms of 
the standard wave functions oc 0„(j/ — ?/o) e"^"" of the 
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nonrelativistic case. For m > the n = eigenmodes 







'My~yo) J ^/2^ 



1 



(A2) 



describe negative-energy electrons with eo_ = —vp m. 
For the parity partner the spectrum and eigenmodes 
are obtained by reversing the sign of m. The n = modes 
are actually the n — 0+ modes Xo+,yo{^) (0; </'o)' which 
describe electrons with positive energy eg , = vp m. 



APPENDIX B: VACUUM ENERGY 

In this appendix we outline the calculation of the vac- 
uum energy in a magnetic field B, discussed somewhere 
below Eq. (|2.26p . Via normal ordering in H, one obtains 
the vacuum energy per unit area 



This is to be compared with the energy density in the or- 
dinary {B = 0) vacuum, e^"^" = (— 5)2wFX]k Vk^ + m"^ 
with the Fermi momentum fcp chosen to give the same 
number of states, = A:|/(27r) = {2N +1)/{2'k£'^), as in 
the B case. The observable vacuum energy density 



thereby becomes finite for N 



N 



oo. 



n=l 



27r^2 



-C(-i)-iVA + ---], 



(B2) 



where A = ^m'^f and -C(-^) ~ 0.2079. The x sector 
also leads to the same Cvac- 



1 

2^ 



1 ^ 
(-^)(2E^" 



\m\vF). 



(Bl) 



n=l 



^ K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and 
A. A. Firsov, Nature (London) 438, 197 (2005). 

^ Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim, Nature 
(London) 438, 201 (2005). 

3 Y. Zhang, Z. Jiang, J.P. Small, M.S. Purewal, Y.-W. Tan, 
M. Fazlollahi, J.D. Chudow, J. A. Jaszczak, H. L. Stormer, 
and P. Kim, Phys. Rev. Lett. 96, 136806 (2006). 

* Y. Zheng and T. Ando, Phys. Rev. B 65, 245420 (2002). 
^ V. P. Gusynin and S. G. Sharapov, Phys. Rev. Lett. 95, 

146801 (2005). 

^ N. M. R. Peres, F. Guinea, and A. H. CastroNeto, Phys. 
Rev. B 73, 125411 (2006). 

K. Nomura and A. H. MacDonald, Phys. Rev. Lett. 96, 
256602 (2006). 

^ J. Alicea and M. P. A. Fisher, Phys. Rev. B 74, 075422 

(2006) . 

® V. P. Gusynin, V. A. Miransky, S. G. Sharapov, and L A. 

Shovkovy, Phys. Rev. B 74, 195429 (2006); see also, E. 

V. Gorbar, V. P. Gusynin, V. A. Miransky, and I. A. 

Shovkovy, Phys. Rev. B 66, 045108 (2002). 
" J.-N. Fuchs and P. Lederer, Phys. Rev. Lett. 98, 016803 

(2007) . 

^ A. N. Redlich, Phys. Rev. Lett. 52, 18 (1984); R. Jackiw, 

Phys. Rev. D 29, 2375 (1984). 
^ A. J. Niemi and G. W. Semenoff, Phys. Rev. Lett. 51, 2077 

(1983). 

^ G. W. Semenoff, Phys. Rev. Lett. 53, 2449 (1984). 

S. M. Girvin and T. Jach, Phys. Rev. B 29, 5617 (1984). 
^ K. Shizuya, Phys. Rev. B 45, 11 143 (1992). 

S. M. Girvin, A. H. MacDonald, P. M. Platzman, Phys. 

Rev. B 33, 2481 (1986). 

M. V. Berry, Proc. R. Soc. London, Ser. A 392, 45 (1984). 

* As in the Tomonaga-Luttinger model, the infinite depth of 



the Dirac sea is an artifact of the continuum model (2.1) 
whose validity is limited to a (far) smaller energy range 
than the original lattice model. It is this energy range that 
sets the scale of the spectrum cutoff, and regularization 
provides a controlled way to derive the (low-energy) re- 
sponse independent of the cutoff. 

N. Fumita and K. Shizuya, Phys. Rev. D 49, 4277 (1994). 
To be precise, it is the real part of the response (such as 
K'^' that enjoys this property; obviously, it is not shared 
by the imaginary part. 

2^ F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 

^■^ To be precise, tj-n becomes insensitive to the transitions 
across the Dirac sea when the four possible types of tran- 
sitions (-A A ± 1) and {~N N ±2) are combined. 
E. Fradkin and F. A. Schaposnik, Phys. Lett. B 338, 253 
(1995); F. A. Schaposnik, Phys. Lett. B 356, 39 (1995). 

D. -H. Lee and S.-C. Zhang, Phys. Rev. Lett. 66, 1220 
(1991). 

2^ K. Shizuya, Phys. Rev. B 63, 245301 (2001); ibtd. 65, 

205324 (2002); Int. J. Mod. Phys. B 17, 5875 (2003). 
2« S. C. Zhang, Int. J. Mod. Phys. B 6, 25 (1992). 
^"^ J. K. Jain, Phys. Rev. Lett. 63, 199 (1989); 

A. Lopez and E. Fradkin, Phys. Rev. B 44, 5246 (1991); 

B. I. Halperin, P. A. Lee and N. Read, Phys. Rev. B 47, 
7312 (1993). 

S. Deser, R. Jackiw, and S. Templeton, Ann. Phys. (N.Y.) 
140, 372 (1982). 
2^ T. Ando, J. Phys. Soc, Jpn. 75, 074716 ( 2006); 

E. H. Hwang and S. Das Sarma, cond-mat/0610561 

B. Wunsch, T. Stauber, F. Sols, and F. Guinea, New. J. 
Phys. 8, 318 (2006). 
^° C. Kallin and B. I. Halperin, Phys. Rev. B 30, 5655 (1984). 



